RG/Pade Estimate of the Three-Loop Contribution to 
the QCD Static Potential Function 

F.A. Chishtie* and V. Eliafl 
Theory Group, KEK, Tsukuba, Ibaraki 305-0801, Japan 



O 

o 

(N 

o 

o 



> 

(N 

in 

o 

o 






X 



Abstract 

The three renormahzation-group-accessible three-loop coefficients 
of powers of logarithms within the MS series momentum-space for 
the QCD static potential are calculated and compared to values ob- 
tained via asymptotic Pade-approximant methods. The leading and 
next-to-leading logarithmic coefficients are both found to be in ex- 
act agreement with their asymptotic Pade-predictions. The predicted 
value for the third RG-accessible coefficient is found to be within 7% 
relative |error| of its true value for n/ < 6, and is shown to be in exact 
agreement with its true value in the nj — > oo limit. Asymptotic Pade 
estimates are also obtained for the remaining (RG-inaccessible) three- 
loop coefficient. Comparison is also made with recent estimates of 
the three-loop contribution to the configuration-space static-potential 
function. 

The perturbative portion of the QCD static potential is presently known to 
two subleading orders of perturbation theory ||l], |^, ^. This potential may be 
expressed as an integral over an MS perturbative-QCD series in momentum 
space. 
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where 



x(/i) = as(/i)/vr, L{n, g ^) = hgin^/q^), 



and where the momentum-space series within the integrand of 
form 



(1) 

(2) 
is of the 



W[x, L] = x[l + (ao + aiL)x + {bo + hL + h2L^)x^ 
+ (co + ciL + 021"^ + c:iL^)x^ + ...] 



(3) 
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with the following known coefficients [Q]: 

ao = 31/12 - 5n//18, ai = 11/4 - n//6, (4a) 

bo = 28.5468 - 4.14714n/ + 25nJ/324, (4b) 

6i = 247/12 - 229n//72 + 5nJ/54, (4c) 

62 = 121/16 -lln//12 + nJ/36, (4d) 

The three-loop order momentum-space coefficients c^ have not been calcu- 
lated. [] It should be noted, however, that the series' convergence may 
be problematical for values of x near 0.1; e.g., if ^/ = 3 and /x^ = q"^, 
W[x,0] = x{l + 1.75x + 16.80x^ + ...]. There is clearly phenomenological 
value in having some knowledge of the next-order coefficients c^ within the 
series (3), even if one chooses sufficiently large values of fi to ensure that the 
expansion parameter x{n) remains small. 

The all-orders momentum space static potential should ultimately be in- 
dependent of the MS renormalization parameter fi: 

^^'^,W[xi^^),Li^^,q')] = 0. (5) 

Eq. (5) corresponds to the following renormalizat ion-group (RG) equation 
for the series iy[x, L]: 



with § 



J 00 

(3(x)^f^'—xifi) = -Y.P,x'+\ (7) 

a/i A:=0 

Po = 11/4 - nf/6, (8a) 

p^ = 51/8 - 19r2//24, (8b) 

P2 = 2857/128 - 5033n//1152 + 325nJ/3456. (8c) 



"'^A leading-log three-loop contribution in configuration space is extracted in refs. [4] 



To leading and next-to-leading orders in perturbation theory, Eq. (6) is 
manifestly satisfied by the known coefficients (4) for the perturbative series 
(3): 

= [-^ + (^^)w[x,L] = {a,-(3,)x' 

+ {hi - 2ao/5o -(3i)x^ + (2^2 - 2(3Qai)x^L 

+ 0{x\x^L,x^L'). (9) 

The coefficients of x^, x^ and x^L in (9) are all seen to vanish for known series 
coefficients (4) and /9-function coefficients (8): i.e., the known values of hi 
and 62 are seen to uphold the RG-equation (6) by satisfying its perturbative 
formulation (9). However, it is important to note that (9) may also be utilized 
to extract all but one of the three- loop coefficients c^ in the series (3). The 
coefficients of x^.,x^L and x^LP' in (9) respectively vanish provided 

ci = 290.769 - 60.4881^/ + 3.2440^5 - 25nJ/648, (10) 

C2 = 1639/16 - 4129n//192 + 377nJ/288 - 5nJ/216, (11) 

C3 = (11/4)^ - 121n//32 + llnJ/48 - nJ/216. (12) 

Consequently, the only RG-inaccessible three-loop-order term in the series 
(3) is cq. This coefficient can be obtained only via a direct perturbative 
calculation, which has not yet been performed. 

In the absence of such a three-loop calculation, we employ the anticipated 
error of Fade approximants in predicting next-order terms of a field theoret- 
ical series in order to obtain an estimate of all four three-loop coefficients 
Cfc within (3). The predictions for {01,02,03} can then be compared to their 
true values (10-12) to check the validity of the estimation procedure. The 
procedure we describe below has already been employed in a large number of 
applications: QCD (3- and 7-functions [0, 0, §, the SQCD /3-function [^, |[, 
QCD current-correlation functions 0, |10|, the renormalizat ion-group func- 
tions of 0(N)-symmetric massive scalar field theory ||, H, |11|, Higgs decays 
1^, |10|, |12|, Higgs-mediated scattering processes |l^, and QCD corrections 



to inclusive semileptonic B-decays |]T^. The general method we employ is 
described in refs. 0, [§ and |T^; we restate its development here for conve- 
nience. 



Consider a perturbative series 

W{x) = l + Rix + R2X^ + R^x^ + ... + Rnx^ + ... . (13) 

For many such series, the series sum can be approximated by an [A^|M] Pade- 
approximant, where N and M are respectively the degrees of numerator and 
denominator polynomials within the approximant. If only the next-to- leading 
term Ri is known, for example, the Pade approximant 

Ty[Q|i](x) = -^ =l + Rix + Rlx^ + ... (14) 

would predict a value of R\ for the coefficient R2 in (13). Somewhat more 
realistically, if Ri and R2 in (13) are both known, the Pade approximant 

W^'\'\x) = ^\^^\~^''^'^'' = l + R,x + R2x' + {Rl/R,)x' + ... (15) 
1 - [R2/ Ri)x 

leads to the predicted value R2/R1 for the coefficient R^ in (13). Generally, 
one finds that the higher the degree of the approximant, the more accurate 
the prediction of the next unknown coefficient of the series will be. Sup- 
pose one now utilises an [A^ — 1 1 1] approximant to estimate the coefficient 
Rn+1 within (13), based upon knowledge of all previous series coefficients 
{i?i, i?2, ••••, -Rat}- For perturbative field-theoretical series, it is often found 
that the relative error in such an estimate is inversely proportional to N 

{r'Z'i - R%Ti) /RZ\ = -A/N. (16) 

The constant A in (16) can be estimated by comparing the [0|1] -approximant 
estimate for R2 {i.e., -Rf'^'^ ' = Rf) against R2^s true value. One then finds 
via (16) that 

A = l- RI/R2. (17) 

In the series (13), let us suppose we only know the subleading and NNLO 
coefficients Ri and R2, as is the case for the series (3). If the [1|1] approximant 
prediction R^^"^ ' = R^/Ri has a relative error described by (16), we can 
substitute (17) into (16) to obtain the "true" value for R^ algebraically M, ITl 



^^ ~ l-A/2~ Rl + R,R2 ^ ^ 



Of course, the validity of this result can be ascertained only by seeing 
how well it predicts coefficients that can be extracted by other means. For 
the case of the series (3), we identify the known coefficients Ri and R2 as 
polynomials in the logarithm L [ffl]: 



R2 



Ri 

bo + biL H 



= ao 4 
b.L' 



aiL = flo + PqL, 
= bo + {2ao(3o + l3i)L + (3lL\ 



(19) 
(20) 



Substituting (19) and (20) into (18), we obtain the following "large-L" series 
expansion for R3,: 



R^ 



where 



+ 
+ 



(ag/3o/2 + 5/3o&o/2 + 5ao/3i/2 + 7PI/AI3o)L' 
(3l{2D,D2 -Dl- D;) + 3/3o(2ao/3o + /5i)(/^? 



D, 



- 3(2ao/3o + PifD,/Po - SPoboDi 

+ (2ao/?o + Pi)VP! + 6&o(2ao/5o + /5i)//?o 



(21) 



Di = {Qao(3o + (3i)/2(3l (22a) 

D2 = [5alPo + boPo + aoPi]/2Pl (22b) 

D; = aoial + bo)/2l3l (22c) 

As is evident from (3), R^ should be a degree-3 polynomial in the logarithm 
L. A direct comparison of equivalent powers of L in (3) and in (21) leads to 
the following predictions: 

cf^'- = (31 (23) 

c^^^'^-=3ao/?o' + 5/5o/?i/2, (24) 

(fr""- = al(3o/2 + 5/5o6o/2 + 5/?iao/2 + 7(3l/A(3o, (25) 

in addition to the predicted equivalence of the unknown coefficient cq with 
the lengthy square-bracketed term in (21). 

^ The formula (18), for example, is surprisingly accurate in predicting the known four- 
loop order /3- function coefficient in 0(N)-symmetric massive scalar field theory [|l3| . 

•^Estimation of higher-order terms via such a series expansion is denoted in ref. [10] as 
the "APAP"' procedure. 



The prediction (23) is in exact agreement with (12), the RG-determination 
of C3, as is evident by substituting (8a) into (23). Surprisingly, the predicted 
value (24) for C2 is also in exact agreement with the RG value (11), as is 
evident from direct substitution of (8a), (8b) and (4a) into (24). Note that 
this agreement for both coefficients is true for all values of Uf, indicating 
that the asymptotic error formula (16) replicates the RG-invariance of the 
series (3) to leading and next-to-leading order in the logarithm L, a most 
surprising result. 

The formula (16) cannot, of course, replicate RG invariance to all orders 
in L, since the infinite series (21) which follows from it is not a degree-3 
polynomial in L. Nevertheless, the coefficient of L in (21) is strikingly close 
to the corresponding coefficient Ci within (3), as obtained via RG- methods 
in (10). In Table 1, such RG determinations of ci are compared to the 
prediction (25). As is evident from the Table, the predicted values for ci 
underestimate corresponding RG values by less than 7% for nj < 6, with 
the best agreement seen curiously to occur at Uf = 6. This feature may be 
understood by noting that the large-n/ behaviour of the estimate (25), 



ci — > -25n^./648, (26) 



is in exact agreement with that of (12), the RG-determination of Ci. 

Table 1 also presents estimates of the coefficient Cq, as obtained from 
the (square-bracketed) L° term in (21). This coefficient, as noted earlier, 
cannot be extracted from lower-order terms via RG-methods. It is never- 
theless encouraging to note that corresponding predictions for C3 and C2 are 
exact, and that predictions of Ci are nearly so. Thus, we have obtained 
in Table 1 asymptotic Pade-approximant estimates for the three-loop coef- 
ficient Co, which, in conjunction with explicit i?G-determinations (10-12) of 
the other three-loop coefficients {ci, C2, C3} occuring within the perturbative 
series W[x, L] (3), constitute a prediction for the full three- loop contribution 
to the static-potential integrand (1). 

The coordinate-space potential corresponding to the series (3) can be 
obtained via (1) through use of the following identities 0: 



(27r)3 p Aur'' 



d^q _gig.r M/^Vi!l _ 2 (log(/xr) + -fE) ^27b) 






(27r)3 g2 4^^ 



A_gig.r[log(yuV^^)] 



(27r)3 g2 

8 (log(/ir) + 7e)' + 27r^ (log(/ir) + 7^) + 16C(3) 

47rr 



(27d) 



where 7^5 = 0.577216 and C(3) = 1.202057. Following ref [16], we set // = 1/r 
and find that 

n (l/r) °° 
Kertir) = ^^^^ E K,<(l/r), (28) 

'^ n=0 

Vo = -4/3, (29a) 

Fi = -(4/37r)[ao + 27£;ai], (29b) 

V2 = -(4/37r2)[6o + 27^61 + (47I + y)62], (29c) 



2 
F3 = -(4/37r=^)[co+27sCi + (47|+y)c2 + (87| + 2n^^E + 16C(3)) C3]. (29d) 

For arbitrary Uf, values of {oq, cti, &0) ^i? ^2} are given by (4), and values 
of {ci,C2,C3} are given by (10,11,12). In Table 2 we display values of the 
coefficients Vi_3 obtained via (29) for Uf = {3,4,5}. The estimate for V3 
is obtained through use of (10, 11, 12) and the estimated values for cq in 
the final column of Table 1. In Table 2 we also list values of V3 estimated 



via renormalon-matching (RM) considerations |jT6[, as well as corresponding 
large-/3o estimates of V3 [[l^ . Striking agreement of all three estimation pro- 
cedures is clearly evident in Table 2. However, it must be noted that V3 is 
not very sensitive to Cq [the only RG-inaccessible coefficient in (29d)] when 



/i = 1/r. If one uses (29d) to extract Cq from V3, for example, one finds that 
the V3 values -38.4, -37.34 (the RM value), and -34.06 (large Pq) tabulated 
in Table 2 for nj = 3 respectively correspond to cq values of 142 (our Table 
1 RG/Pade estimate), 116, and 40. 

An alternative approach to estimating cq follows from a least-squares 
fit of the asymptotic Pade-approximant prediction (18) to the three-loop 
momentum-space contribution's explicit dependence on L, 

R3 = co + ciL + csL^ + c-sL^, (30) 

over the entire ultraviolet (/x^ > g^) region, a procedure which has been 
employed previously in a number of different applications []T2], |T3], |T^]. If we 



define w = q'^/fJ.'^ [i.e., log{w) = —L], such a procedure entails optimization 
of 



Jo 



2 H'^diA 

Co + ci log(w;) - C2 log^(w) + C3 log^(w) 



Rf{w) + Ri{w)R2{w] 

(31) 
with respect to Cq, where 

R^{w) = ao-ai \og{w), i?2(w) = 60-^1 log(w) + ^2 log^(w), (32) 

and where the set of known coefficients {oq, ai, 6o; ^i; ^2, Ci, C2, C3} is given by 
(4) and (10-12). Unlike previous applications in which large-L expansions of 
(18) are quite consistent with least-squares fits, such a fit is seen to lead to 
values of cq that are ~50% larger than those of Table 1: 

nf = 3: x^N =42679-405.8co + cg ^co =203, (33) 

X^N = 22142 - 291. Ico + eg ^ cq = 146, (34) 

;^2[co] = 9501 - 189.9co + eg -^ Cq = 95, (35) 

x'[co] = 2901 - 104.7C0 + cl ^ Cq = 52, (36) 

In assessing the accuracy of (33-36), it should be noted that such least- 
squares fitting could also be employed to fit simultaneously a// /ottr three-loop 



'/ = 


4 : 


% = 


5 


^/ = 


6 



For example, in semileptonic b —> u decay, the n/ = 4 cq coefficient has a large- L- 
expansion value of 166, in approximate agreement with the estimate Cq = 188 obtained 
in [Q via (28) with known values ||l^ for {ao, ai, 60, &i, fe} and RG-determinations [Q 
of {ci, 02,03} appropriate for b — > u£~'i>e- 



coefficients {cq, Ci, C2, C3}, as has been done before in a number of applications 
[1T2| , p!3| , in which the fitted values for {01,02,03} closely approximated their 
known RG values. However, such a procedure completely fails for the series 
(3): when nf = 3, optimization of (31) with respect to cq-s yields values 
[co = 258, ci = 54.7, C2 = 66.3, C3 = 10.2] that differ substantially from true 
values [ci = 137.46, C2 = 49.078, C3 = 11.391] obtained from eqs. (10,11,12). 
A similarly large estimate of Cq for the nf = 3 case is obtained directly via 
(18) in the L = (small-log) limit, in which R2 = &o (4b) and Ri = oq 
(4a). Such an approach yields cq = 273, a value quite comparable to that 
obtained above (cq = 258) by simultaneous least-squares fitting of all four 
three loop coefficients Co_3. If one inputs this small-log estimate for cq (= 273) 
into (31) and then minimizes with respect to the RG-accessible coefficients 
Ci_3, the estimated values for these coefficients will be even worse than those 
characterising the full least squares fit (ci ~ 33, C2 — 74, C3 ~ 9.6) - values 
which are inconsistent (particularly ci) with the RG determinations of these 
same parameters. 

Such discrepancies suggest that the large-L {i.e. large fi or short distance) 
Co estimates of Table 1, which reproduce exact RG values for C2 and C3 and 
closely approximate RG values for ci, be taken more seriously than either the 
Co estimates (33-36) obtained via least-squares fitting over a broad range of //, 
or other {e.g. small-log) approaches to estimating Co within an asymptotic 
Pade-approximant context. It is evident that such alternative asymptotic 
Pade-approximant estimation procedures are of little value if not tied to some 
way of successfully estimating ci_3, the RG-accessible three loop coefficients. 
By this criterion, the large-L estimates of Table 1 have the most substantial 
credibility. 

However, it is important to remain cognisant of the relative insensitivity 
(noted earlier) of the configuration-space static potential to the parameter 
Co at its benchmark /i = 1/r length scale. At this scale, RG-accessible co- 
efficients alone dictate that the n/ = 3 three-loop contribution V3 is given 
by 

V3 = -4 [co + 752.0] /Sn^ (37) 

a result which follows from substitution of (10-12) into (29d). Since all 
nf = 3 estimates of co, as delineated in the previous two paragraphs, are small 
compared to 752.0, V3 is surprisingly insensitive to this unknown parameter; 
the four disparate estimates 142 [large-L], 203 [eq. (33)], 258 [full least- 
squares treatment], and 273 [small-log] all correspond to V3 values near -40 



[-38.4, -41.1, -43.4, and -44.1, respectively]. Thus, even if one trusts Pade 
methods only to give a factor-of-two-accuracy estimate for the magnitude 
of Co, RG-accessible coefficients alone are sufficient to extract a surprisingly 
concise range for the three- loop contribution V3. In other words, the n = 
1/r estimate for V3, the three-loop quantity ultimately of phenomenological 
interest to us, is subject to substantially less theoretical uncertainty than the 
parameter cq. 

A final and necessary caveat, however, is the possibility that new diagram- 
matic topologies (and their corresponding group theoretical factors) known 
to enter the QCD static potential at three-loop order |^ may further cir- 
cumscribe the applicability of using lower-order terms to predict the three- 
loop contribution cq, as in a Pade-approximant approach. This situation 
is entirely analogous to the theoretically uncertain light-by-light scattering 
contributions known to enter the muon's anomalous magnetic moment at 
sufficiently high order, as well as the quartic Casimir terms first appearing in 
the QCD /3-function series at four-loop order. Pade-approximant based tech- 
niques cannot be expected to predict terms characterised by new higher-order 
group-theoretical factors [0. Nevertheless, such contributions do not neces- 
sarily dominate the first order in which they appear, nor do they necessarily 
devalidate Pade estimates for that order. For example, the asymptotic Pade- 
approximant estimate of the {N^ = 3) four-loop contribution to the QCD 
/3-function /?3 = 23600- 640lMf+ 350nj+ 1.4993171) (estimated numbers 
are italicised) is in quite reasonable agreement with the exact result [|T^] /^s = 
29243.0 - 6946.30n/ + 405.089n^ + 1.49931ni^, though in much closer agree- 
ment with the calculated result with quartic Casimir terms excised : /^s = 
24633 - 6375nf + 398. Sn^ + 1.49931raj. Thus, based on the limited informa- 
tion available, the remarkable success of the asymptotic Pade-approximate 
large-L expansion in predicting those three-loop momentum-space static po- 
tential terms that are also extractable by RG methods encourages some con- 
fidence in the corresponding large-L prediction of the RG-inaccessible pa- 
rameter Co (modulo the above-mentioned uncertainties characterising Pade 
approaches), at least for purposes of predicting V3, the three-loop contribu- 
tion to the configuration- space static potential. 
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Pred 


\^^Pred_^ROy^m 


Pred 

Co 





290.77 


272 


6.3% 


313 


1 


233.49 


218 


6.5% 


250 


2 


182.46 


170 


6.8% 


193 


3 


137.46 


128 


7.0% 


142 


4 


98.251 


91.4 


6.9% 


97.5 


5 


64.606 


60.6 


6.2% 


60.1 


6 


36.291 


35.3 


2.8% 


30.5 



Table 1: Comparison of predicted and RG values for the three-loop coefficient 
Ci. Also displayed are predicted values for the RG- inaccessible coefficient Cq. 
Note that these cq estimates are the same sign and approximate magnitude 
as the RG values for ci listed in the second column. 



^f 


Vi 


V2 


V3 


V,^'^ 


yWo 


3 


-1.84512 


-7.28304 


-38.4 


-37.34 


-34.06 


4 


-1.64557 


-5.94978 


-28.7 


-27.63 


-27.03 


5 


-1.44602 


-4.70095 


-20.5 


-19.46 


-21.05 



Table 2: Configuration-space coefficients (28) of the configuration-space 
static potential. The column labeled V3 is obtained using RG-determinations 
of Ci, C2, C3 and the Table 1 estimate of cq within eq. (29d). The column la- 
beled Vg^^^ is obtained from eq. (22) of ref 
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The column labeled K 



lists large- /3o estimates [0 that are also tabulated in ref. |jT6| 
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